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Abstract. The merger dynamics of a black hole-neutron star (BHNS) binary is influenced 
by the neutron star equation of state (EoS) through the latter' s effect on the neutron star's 
radius and on the character of the mass transfer onto the black hole. We study these 
effects by simulating a number of BHNS binaries in full general relativity using a mixed 
pseudospectral/finite difference code. We consider several models of the neutron star matter 
EoS, including T = 2 and V = 2.75 polytropes and the nuclear-theory based Shen EoS. 
For models using the Shen EoS, we consider two limits for the evolution of the composition: 
source-free advection and instantaneous /3-equilibrium. To focus on EoS effects, we fix the 
mass ratio to 3:1 and the initial aligned black hole spin to a/m = 0.5 for all models. We 
confirm earlier studies which found that more compact stars create a stronger gravitational 
wave signal but a smaller postmerger accretion disk. We also vary the EoS while holding the 
compaction fixed. All mergers are qualitatively similar, but we find signatures of the EoS in 
the waveform and in the tail and disk structures. 



PACS numbers: 04.25.dk, 04.40.Dg, 04.30.Db, 47.75.+f, 95.30.Sf 

1. Introduction 

Black hole-neutron star (BHNS) inspirals and mergers are excellent gravitational wave 
sources. Also, a BHNS merger may leave a hot, massive accretion disk around the black hole 
(BH), a promising setup for producing a short-duration gamma-ray burst (GRB). However, 
reliable predictions regarding the waveform and post-merger state of any specific BHNS 
merger can only be obtained from fully relativistic numerical simulations. 

The merger is strongly affected by the NS equation of state (EoS), i.e., the dependence 
of the fluid pressure on density, temperature, and composition. The EoS determines the radius 
i?NS ofaNS of given mass Mns, and thus it also fixes the star's compaction C = Mns/-Rns- A 
larger star will be tidally disrupted by the BH's gravity at a larger distance. Also, because the 
EoS affects dR^s/dM^s, it can influence the character of the mass transfer from the star to the 
BH, with mass transfer tending to be more stable for stiffer EoS. Differences in the stability of 
mass transfer could conceivably result in qualitative differences in the merger dynamics. The 
NS EoS at low densities is known, but at high densities it is not well constrained by theory, 
experiment, or observation. It must, therefore, be treated as another parameter in BHNS 
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simulations. If the EoS significantly affects the waveform, then comparison of numerically 
generated waves with observations may someday provide information about the EoS. 

Attempts have been made to estimate the effects of NS EoS in the context of Newtonian 
BHNS simulations. Newtonian simulations have considered poly tropic EIU[3lH, Lattimer- 
Swesty (51 [6l and Shen EKHHH EoS. These simulations showed large qualitative differences 
for different EoS assumptions. For Lattimer-Swesty nuclear matter, the NS disrupts in one 
mass transfer event, and a large post-merger disk is created. For Shen nuclear matter, a NS 
core can survive multiple mass transfer events, and the postmerger disk is much smaller. 
There are indications that the differences are not quite so dramatic when general relativity 
(GR) is accounted for. The use of GR-mimicking potentials |[T0l [TT| tends to eliminate 
episodic mass transfer lfT2l[T3l . Simulations of large mass-ratio cases using the conformally 
flat approximation of GR also found mass transfer to be less stable (and, therefore, surviving 
cores to be less likely) than in Newtonian evolutions |fl4|. These conformally-flat-GR studies 
only considered fairly soft EoS, however. 

Published BHNS simulations in full GR have thus-far restricted themselves to modeling 
NS as T = 2 polytropes (TH EES El M\ EES IM EQ. These have investigated the effects of 
varying one important EoS-related variable, namely R^s flTVl UHl [2TJ. Varying R NS between 
12 and 15km for a fixed M B h and M NS , they find, unsurprisingly, that larger stars disrupt 
farther from the BH, leading to larger disks but weaker merger waveforms. These simulations 
could not address the importance of the stiffness of the EoS in affecting the mass transfer. 

Another important question is whether BHNS mergers eject significant amounts of 
NS matter (perhaps including r-process elements (22l 12310 into the interstellar medium. 
Newtonian (6]|, pseudo-Newtonian lfT2l[T3l . and conformal-GR lfi~4l simulations predict large 
ejecta masses, but this has not yet been seen in full GR. However, microphysics in the tidal 
tail could have important effects. As the density decreases, nucleons recombine into nuclei, 
heating the gas through the released binding energy. The resulting increase in thermal pressure 
can strongly affect the tail Q, perhaps helping to unbind material. Material may also be 
ejected in a neutrino-driven or magnetic field-driven disk wind, effects which are also not 
captured in current GR simulations. 

In this paper, we study the imprint of the NS EoS on the merger dynamics and the 
resulting waveform and disk. Our investigation is based on simulations of BHNS mergers 
in full GR carried out with the Cornell-Caltech code SpEC. In a previous paper [19], we 
demonstrated this code's ability to evolve BHNS inspirals and mergers. Since then, we have 
made a number of improvements which reduce the constraint violation during merger by an 
order of magnitude. For this study, we fix the binary mass ratio to 3:1 and the initial BH 



spin to be 



Si 



BH 



/M B h 2 = 0.5 (aligned with the orbital angular momentum), and we vary 
the assumed EoS. We evolve n — 1 polytropes using a T = 2 Gamma-law and n = 4/7 
polytropes using T = 2.75. We also perform simulations using the tabulated, nuclear-theory 
based Shen EoS (71 HJ. We vary the EoS and the compaction separately to isolate the effect of 
each. 

All our mergers are qualitatively similar: we find no cases of episodic mass transfer, 
no instances of measurable unbounded outflow, and no cases in which the disruption fails 
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to produce a significant disk. There are, however, some quantitative EoS signatures in the 
waveform above 1kHz and in the disk. We confirm earlier findings regarding the effects of 
NS compaction. For a given compaction, stiffer EoS produce larger, longer-lived tidal tails. 
In every case, the post-merger disk has a mass of 0.05-0. 1M Q and an average temperature of 
order an MeV. 

In Section 2, we describe important improvements to our code from |fT9l . In Section 3, 
we describe the models we evolve. In Section 4, we present the results of our evolutions. In 
Section 5, we summarize our findings and consider what remains to be done to accurately 
sample all of the interesting regions of BHNS parameter space. 

2. Numerical methods 

2.1. Code improvements 

The main improvements to our code from |fT9ll come from a better allocation of grid points 
when solving the fluid equations, modified gauge conditions and atmosphere prescriptions. 

Our simulations use two grids |fT9ll : a pseudospectral grid on which Einstein's equations 
are solved, and a finite difference (FD) for the relativistic fluid equations. For efficient 
evolutions, the FD grid should adapt to the configuration of the fluid. One approach is to 
use adaptive mesh refinement, as in some FD-based GR codes E4l l25l l26l . We choose 
instead to modify the mapping between the two grids as soon as a significant flow of matter 
is approaching the boundary of the FD grid. A translation and rescaling of the coordinates 
allows the grid to move, grow, and shrink so as to follow the fluid evolution. To increase the 
resolution in the neighborhood of the black hole, we use a map of the form 

r' — r (r < Rq) 

r' = ar 3 + br 2 + cr + d (R < r < (1) 
r' = ar (r > Ri), 

where r is the coordinate distance to the center of the hole, R and Ri are predetermined 
length scales, the parameters (a, b, c, d) are chosen so that the map is C 1 , and a varies with 
the size of the FD grid. (Once the grid becomes large, we maintain a constant resolution in 
the region r < R and vary a to fix the location of the outer boundaries.) 

The gauge in the generalized harmonic formulation is set by specifying the functions 
H a = g bc T abc . During the inspiral, we fix H a in the moving frame, as in our earlier paper |[T9ll . 
During the merger, we continue to hold H a fixed in the moving frame in the region near the 
black hole, and we exponentially damp H a to zero in the region far from the hole. We find 
that fixing H a near the excision zone significantly reduces constraint violations. 

As Faber et al. 11271 have pointed out, the inversion from conservative to primitive 
hydrodynamic variables is only possible if S = g tj SiSj < S mSLK = r(r + 2D). After each 
evolution step, we impose the condition S < = fS maiX . We find 0.99 < / < 1 

is necessary to avoid causing large effects on the evolution of the tidal tail. Once the 
values are "fixed" in this way, primitive variables can be reconstructed, but they may still 
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be unreasonable in very low density "atmosphere" regions. Therefore, we next apply limits 
on the conformal 3-velocity itj and the temperature T. We emphasize that these limits are 
only applied to low density regions, several orders of magnitude sparser than the star, the 
disk, or the tidal tail. We have checked that our evolutions are insensitive to variations in 
these atmosphere ceilings. 

With our standard resolution, the normalized constraint violations peak at ~ 1% during 
mergers. Convergence tests on the Y = 2.75 case indicate errors in our reported disk masses 
of < 10%. 



2.2. Use of tabulated p/T/Y e -dependent EoS 

For composition-dependent EoS, there is a new independent variable to be evolved: the 
electron fraction Y e . Its evolution equation in conservative form is 

d t (DY e )+d i (DY e v i ) = S v , (2) 

where S u is the source term set by weak interactions and neutrino radiation, effects not 
modeled in our code. Here, we consider two limiting cases. First, we assume the weak 
interaction timescales are much longer than the merger timescale. Then we may set S u = 
and evolve a continuity equation for DY e . As another limit, we assume that weak interactions 
act sufficiently quickly to instantaneously enforce /^-equilibrium. Thus, for a given density p 
and T, Y e is set to the value that makes fi n = fi p + /i e , where jj x is the chemical potential of 
particle x. (We assume [i v is negligible.) This effectively removes Y e as a dynamical variable. 
In neither case do we account for energy loss by neutrino emission. Given that the cooling 
timescale of the disk will probably not be less than about 0.1 second (see, e.g. 11281 ). and our 
simulations last ~10ms, ignoring neutrino cooling is reasonable. 

We use a tabulated EoS with baryon component taken from Shen et al [0 and with 
lepton and photon contributions added. (See ll29l for details.) For low T (< lOMeV), the NS 
EoS is fairly soft (r « 4/3) at low p (< 10 12 g cm" 3 ) and stiffer at high p (T ps 2.75 in the 
Shen model). The Shen EoS predicts a TOV maximum gravitational mass of 2.2M . 



3. Cases 

In Table [Q we present the initial data for the cases we evolve. For this study, we do not 
consider the effects of the binary mass ratio q = Mbh/Msts and set q = 3 throughout. For 
each case, we set the NS baryonic mass to 1.55M . We wish to study cases that lead to 
massive disks, and so we include in each case an initial BH spin of s = Sbu A^bh = 0.5 
orthogonal to the orbital plane. Finally, we neglect the NS spin altogether ana focus solely 
on what is thought to be a good approximation for the most likely scenarios: irrotational 
stars II301I3TTI. 

To study the effect of compaction, we ran two cases with T = 2.75, one with C = 
M NS /R NS = 0.146 and one with C = 0.173. These runs are labeled 'T2.75c.15" and 
'T2.75c.17" in the tables and figures below. For our chosen rest mass, the NS gravitational 
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Table 1. Initial data for all runs. d; n it is the initial coordinate separation and p\ A is the central 
density of the star in units of 10 14 g cm" 3 assuming a star of baryonic mass M^ s = 1.55M Q . 
e = -^j is the eccentricity and B is obtained by fitting the evolution of the trajectory to 



d = A + 


Ait + Bsm(ut + 4> 


). M is the ADM mass of the system. 




Case 


EoS 




q 


Mn.it 




dtait/M 


Pl4 


e 


T2 


T = 2.00 


0.5 


3 


0.144 


0.041 


7.48 


7.1 


0.006 


T2.75c.15 


T = 2.75 


0.5 


3 


0.146 


0.041 


7.45 


4.5 


0.003 


T2.75c.17 


T = 2.75 


0.5 


3 


0.173 


0.041 


7.42 


7.6 


0.014 


Shen-/? 


Shen 


0.5 


3 


0.147 


0.041 


7.59 


5.0 


0.004 


Shen-Adv 


Shen 


0.5 


3 


0.147 


0.041 


7.59 


5.0 


0.004 



mass is about Mns = 1.4M , so the two compactions correspond to radii of 14.4km and 
12.1km. 

To separate the effects of the EoS from those of the compaction, we use three equations 
of state with the same -Rns (C = 0.15). This comparison has not been attempted in previously 
published full-GR numerical studies. We use two T-law EoS: T = 2 (run 'T2") and T = 2.75. 
We also use the tabulated Shen EoS. As mentioned above, we evolve the initial data with 
Shen EoS in two ways: assuming instantaneous /^-equilibrium (run "Shen-/3") and assuming 
S u = (run "Shen-Adv"). 

All the initial configurations are generated using our multidomain spectral elliptic 
solver ll32l to solve for quasi-equilibrium configurations of BHNS binaries in the extended 
conformal thin sandwich formalism [33J. We do not assume conformal flatness, but instead 
choose a conformal metric approaching Kerr in the neighborhood of the BH Il34ll33l . For the 
polytropic runs, the initial temperature T init is zero. For the Shen runs, T init = O.lMeV, and 
the initial Y e is set by assuming /3-equilibrium. To ensure that the comparison is not affected 
by the initial eccentricity of the binary, we also apply the eccentricity removal technique 
devised by Pfeiffer et al. 11351 until e ~ 0.01. (Without eccentricity removal, we would have 
e = 0.01 - 0.06.) 

4. Results 

4.1. Qualitative features of the mergers 

For each configuration, the binary passes through 2-3 orbits of inspiral before reaching a 
coordinate separation of d ^40km (5M), at which point matter starts flowing from the star 
to the BH. Most of the core of the star is then rapidly accreted onto the hole [within t ~3ms 
(100M)] while a large tidal tail forms. For every EoS, even the stiffest (T = 2.75), the star is 
disrupted in one extended mass transfer event. This confirms earlier indications that episodic 
mass transfer does not happen in GR for the realistic range of neutron star EoS. All tails but 
the one formed by the most compact star extend to distances d >400km. All of this material, 
however, remains bound to the system. This result is significant because our Shen-EoS runs 
do account for recombination effects that occur for p < 10 14 g cm -3 . It would seem that it 
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Table 2. The properties of the BH and accretion torus at time t = i mcrgcr + 8.3ms. Wkick is the 
BH kick velocity in km s" 1 . Mdisk is the baryonic mass outside of the black hole, rdisk is the 
radial extent of the disk, p™{^ 12 is the maximum density in units of 10 12 g cm" 3 , (T)di s k is 

2 2 

the density-weighted average temperature in MeV and W 4 ' is the amplitude of gravitational 
waves in the (2,2) mode extracted at r = 75M. Note that the disk continues to evolve at late 
times as rdisk and (T)di s k slowly increase. 



Case 


Mbh/M 


Sfinal 


^kick 


M disk /M NS 


raisk/ 'M 


„max 
Pdisk,12 


(T)disk 




T2 


0.96 


0.69 


83 


0.08 


12 


1 


1.3 


0.012 


T2.75c.15 


0.94 


0.69 


45 


0.13 


12 


1 


0.6 


0.010 


T2.75c.17 


0.96 


0.70 


75 


0.02 


6 


0.3 


1.2 


0.030 


Shen-/3 


0.97 


0.79 


63 


0.07 


10 


1 


2.7 


0.015 


Shen-Adv 


0.97 


0.80 


61 


0.07 


11 


1 


2.5 


0.016 



is the inclusion of GR, and not inferior microphysics, that causes GR simulations to see no 
ejecta. (It is, however, possible that some very low-density material is ejected and that this is 
suppressed in our simulations by the atmosphere prescription.) 

Material from the tail falling back towards the hole has enough angular momentum to 
avoid being immediately accreted: once the core of the star has fallen into the black hole, the 
accretion rate drops and a disk forms from the remains of the tail. This disk is at first neither 
axisymmetric nor time-independent: the density peaks strongly at the junction of the disk 
and the tidal tail, and the matter keeps expanding away from the hole. As the accretion rate 
decreases, a gap opens between the disk and the black hole. At the end of our simulations, the 
outer disk is still quite nonaxi symmetric, and the tidal tail is settling back onto the disk. The 
inner disk has a density p ~ 10 11 — 10 12 g cm -3 and temperature T ~ 0.1 — 10 MeV. It is 
also relatively thick (i/disk/-Rdisk ~ 0.1 — 0.3). 

4.2. Effects of composition, compaction, and stiffness 

The effects of the EoS on the gravitational wave spectrum and on the post-merger disk 
mass are shown in figures [TJ and [2l The leading contribution from the EoS is related to 
the compaction of the star. As Shibata et al. lfT8l EU found for Y = 2 stars, we observe 
that a higher compaction leads to stronger gravitational waves and a smaller disk. The more 
compact star reaches the ISCO with minimal distortion and falls nearly entirely into the black 
hole, with only a small fraction of its mass sent out in a relatively short tidal tail extending 
to d ~250km. The tail forms a small disk of mass ~ 0.02Mtvs. The cutoff frequency of 
the gravitational waves (/ ~ 2kHz) is higher than for more extended stars (/ ~ 1kHz for 
r2.75c.15). 

By comparison, modifying the EoS while keeping the compaction constant seems to have 
more modest effects. For NS compaction of C = 0.15, every EoS predicts a final disk mass 
~ 0.05 — 0.1M NS , i.e. M dis k ~ 0.1M o . The density-weighted average temperature of the 
disk is l-2MeV, although some regions in the disks reach T « 20MeV. Differences in the 
gravitational wave signal are more visible: the T = 2.75 star has a lower cutoff frequency 
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Figure 1. The effect of EoS on the disk mass. Mq is the total baryonic mass outside the BH. 

and peak amplitude. But even then, the differences remain smaller than those due to the 
compactness of the star. The close similarity of the Shen and T = 2 waveforms is surprising. 
A partial explanation may lie in the fact that the Shen-EoS is actually very soft at the densities 
of the matter after the NS is disrupted. In these regimes, Shen is closer to T = 2 than to 
T = 2.75. 

Another strong effect of the stiffness of the EoS is the different behaviors of the tidal 
tails, which can lead to important effects on the post-merger disk dynamics. For all EoS, our 
simulations end with 0.06 — 0.07Mns of matter within 200km of the BH — either in the disk or 
in the process of joining it. The mass and size of the tidal tail, however, vary: for the T = 2.75 
star, the tail mass is 0.06Mns and more than 2% of the matter will eventually go as far as 
2000km away from the BH before falling back on a timescale of 200ms (assuming geodesic 
motion). When the tail is maximally extended, our grid would be unable to adequately resolve 
the whole system. Therefore, we stop evolving when the tail is still expanding. On the other 
hand, the T = 2 EoS leads to a much smaller tail, with a mass of 0.02M NS . Only a negligible 
amount of matter reaches a distance of 500km, and most of the matter would fall back on the 
disk within 25ms of the disruption. For both polytropes, interactions between the disk and the 
tail are strong enough to keep the disk from settling to an axisymmetric state over the duration 
of our simulation. The Shen case also has a tail of mass 0.02Mns, but it is thicker and falls 
back more rapidly onto the disk. The disk will settle more quickly, and is much closer to 
axisymmetry at the end of our simulation. 

For both Shen-/? and Shen-Adv runs, the density-weighted average of Y e remains close 
to (Y e ) = 0.09 throughout the inspiral, so the inspiral is essentially identical for the two runs. 
The mergers are also similar, although the Shen-/3 waveform decays somewhat more rapidly. 
The composition of the disk, however, is radically different. For Shen-Adv, the final (Y e ) is 
about 0.09. (It decreases slightly, even though there are no Y e source terms, because the more 
highly-leptonized central region of the NS is swallowed by the BH.) The baryon mass is about 
86% free neutrons, 7% heavy nuclei, 5% free protons, and 2% alpha particles. The nuclei have 
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Figure 2. The effective amplitude, as defined in Eq. (41) of 1211 . for each run. The assumed 
distance from the source is lOOMpc. 



an average (A) ~ 95, (Z) « 30. If /3-equilibrium is enforced, (Y e ) increases to 0.2 as the 
NS matter decompresses. The baryon mass is about 65% free neutrons, 25% heavy nuclei 
({A) « 75, (Z) « 30), 9% free protons, and 1% alpha particles. The disks produced by the 
Shen-/? and Shen-Adv mergers have similar densities and temperatures: the density-weighted 
average temperature is about 2.5MeV for Shen-Adv and 2.7MeV for Shen-/?, although in each 
case the maximum temperature reaches rs 12MeV, and the average is slowly increasing. 

4.3. The final black hole and disk state 

The final state of the black hole is given in Tabled The final spin s final is in the range 0.7-0.8, 
with the Shen runs having higher Sfl na i because they merge somewhat more quickly and have 
less time to radiate angular momentum. The BH kick velocity is ~10-100km s" 1 

In figure [3l we plot the profiles of the density p, specific entropy s, and specific angular 
momentum j = u^ju t for run Shen-/? at final time t =10ms after the merger. By 3ms after 
merger begins, a distinct torus forms around the BH. By 5ms after merger, matter begins 
to clear out in the region near the hole, and the accretion rate Mdisk drops to a low value 
(r = M disk /M disk w 24ms). The torus has a maximum density of p = 10 12 g cm" 3 located at 
a coordinate radius of r = r c = 30km from the BH. The width of the torus is w 20km, and 
the height is w 10km. At this time, the disk itself has not yet stabilized. The gas more than 
20km from the BH has nearly constant s and j, while the gas nearer the hole has somewhat 
higher s and much lower j. Subsequent motions in the fluid lead to modest but positive radial 
gradients for s and j in most of the high-density region, the exception being the small but 
negative entropy gradient that persists around r w 30km. In the still- settling, lower-density 
outer regions, s and j drop significantly. The angular velocity Q decreases with r everywhere: 
Q ~ r -1 - 2 in the high-density region. The assumption of ,3-equilibrium is probably not good 
in the outer regions of the disk, which are optically thin to neutrino emission. The profiles of 
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Figure 3. The density p, specific angular momentum j, and specific entropy s of the Shen-/3 
disk as a function of cylindrical radius r, shown 10ms (370A/) after the merger. Each point is 
a p-weighted average over the angular and vertical directions. 



p and j are quite similar for the Shen-Adv disk, so these variables appear to depend weakly 
on Y e . 

The longer-term evolution of the disk depends on physical processes not included in 
these simulations. Since dVt/dr < 0, our disks are subject to the magnetorotational instability 
(MRI) ll36l . Turbulence induced by this instability might have the effect of an a-viscosity of 
magnitude a ~ 0.01 — 0.1 IlTTI . This will heat the disk, redistribute angular momentum, and 
drive accretion. The disk's only significant cooling mechanism is neutrino emission 113811391 . 
3D Newtonian simulations predict that a neutrino-cooled disk with a ~ 0.1 and Mdi S k similar 
to ours will accrete on a timescale of ~0. 1 second and release energy in neutrinos at a rate 
L v ~ 10 53 erg s -1 , possibly making such a merger remnant a viable GRB candidate [40]. 

5. Conclusions 

We have investigated the influence of the equation of state on BHNS binaries. We find that 
the NS compaction has a strong influence on the disk mass and the cutoff frequency of the 
gravitational waveform. The effects of EoS stiffness for a fixed compaction are weaker. The 
overall behavior of the merger was found to be independent of EoS stiffness — we never find 
episodic mass transfer or unbound ejecta. The stiffness of the EoS in the lower-density outer 
layers of the NS does affect the merger dynamics, with stiffer EoS leading to larger and 
longer-lived tidal tails, and this effect manifests itself in the merger waveform. We find that 
the evolution of Y e weakly influences the waveform and the final disk mass. 

Our simulations suffer from two classes of limitations. First, we have considered only 
a very small sample of the interesting BHNS parameter space. It is possible that the effects 
of EoS stiffness are more pronounced at different binary mass ratios or different BH spins. 
Also, we have considered only three EoS. A more systematic approach would be to use a 
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single EoS with adjustable parameters (e.g. II4TI I5T0. An adequate EoS for this purpose would 
have to cover the range of likely NS stiffness and radius while capturing all of the important 
physical EoS features, including (for the post-merger evolution) its complicated temperature 
and composition dependence. Other important areas for improvement involve our treatment of 
the NS microphysics. To evolve the final disks realistically, the effects of neutrino radiation 
on the temperature and composition of the matter must be included. The evolution of the 
disk is also strongly affected by the presence of magnetic fields, and in particular by the MRI. 
Simulations that include radiation and magnetohydrodynamics are needed to assess the ability 
of these disks to produce GRBs. 
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